Applications of density matrices 
in a trapped Bose gas 



K. Ch. Chatzisavvas* S. E. Massen^, Ch. C. Moustakidis* C. P. Panos§ 
Department of Theoretical Physics, 
Aristotle University of Thessaloniki, 
54124 Thessaloniki, Greece 

January, 2006 



Abstract 

An overview of the Bose-Einstein condensation of correlated atoms 
in a trap is presented by examining the effect of interparticle corre- 
lations to one- and two-body properties of the above systems at zero 
temperature in the framework of the lowest order cluster expansion. 
Analytical expressions for the one- and two-body properties of the 
Bose gas are derived using Jastrow-type correlation function. In ad- 
dition numerical calculations of the natural orbitals and natural oc- 
cupation numbers are also carried out. Special effort is devoted for 
the calculation of various quantum information properties including 
Shannon entropy, Onicescu informational energy, Kullback-Leibler rel- 
ative entropy and the recently proposed Jensen-Shannon divergence 
entropy. The above quantities are calculated for the trapped Bose 
gases by comparing the correlated and uncorrelated cases as a function 
of the strength of the short-range correlations. The Gross-Piatevskii 
equation is solved giving the density distributions in position and mo- 
mentum space, which are employed to calculate quantum information 
properties of the Bose gas. 
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1 Introduction 

The first theoretical prediction of the famous phenomenon known as Bose-Einstein 
condensation (BEC) was made in 1924 and 1925 by Bose and Einstein [2], 
respectively. In a system of particles obeying Bose-Einstein statistics where the 
total number of particles is conserved, there should be a temperature below which 
a finite fraction of all the particles condense into the same one-particle state [21111 El 
inilZllHj- Seventy years later, in a remarkable experiment, Anderson et al. [H] have 
cooled magnetically trapped ^"^Rb gas to nanokelvin temperatures, and observed 
the BEC. This discovery has generated a huge amount of theoretical investigations 

[iniiniiniiiiiiiiiiiaiiniiniiiHi. 

The main feature of the trapped alkali-metal and atomic hydrogen systems 
(which obey the Bose-Einstein statistics) is that they are dilute. The crucial 
parameter defining the condition of diluteness is the gas parameter x = '^o.^ > where 
n is the density of the system and a is the s-wave scattering length |14j . There are 
two ways to bring x outside the regime of validity of the mean field description. 
The first one is to increase the density, while the second one to change the effective 
size of the atoms. 

The diluteness of the gas is ensured when the effective atomic size is small 
compared both to the trap size and to the interatomic distance. However, the 
effects of inter-particle interactions are of fundamental importance in the study of 
the BEC dilute-gas where the physics should be dominated by two-body collisions 
described in terms of the s-wave scattering length a. In the case of positive a, 
it is equivalent to consider a very dilute (atomic) system of hard spheres, whose 
diameter coincides with the scattering length itself jl3j . The natural starting point 
for studying the behavior of those systems is the theory of weakly interacting 
bosons which, for inhomogeneous systems, takes the form of the Gross-Pitaevskii 
equation |19l l2Uj . This is a mean- field approach for the order parameter associated 
with the condensate 1^. 

In the present work we study BEC in a phenomenological way where the Bose 
gas is considered as a many-body system j2^. In particular, we study the ground 
state of a system of correlated bosonic atoms at zero temperature, trapped by a 
harmonic oscillator potential (HO). The key quantities for this effort is the one- 
and two-body density matrices j^. As the mean-field approach (non-interacting 
atoms) fails to incorporate the interparticle interactions which are necessary for the 
description of the correlated Bose system, we introduce the repulsive interactions 
among the atoms, through the Jastrow correlation functions / (|ri — r2|) ^23; . 
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We focus our efforts to the calculation of the one- and two-body density and 
momentum distributions and the calculation of the static structure factor [231 
1241 1251 One-body density and momentum distributions are complementary 
descriptions of the Bose gas and related directly with the mean-square radius and 
mean kinetic energy of the trapped Bose gas respectively. In addition the two-body 
density distribution is related with the calculation of the static structure factor, 
a quantity which gives information for the ground and excited states of the gas. 
Special effort has been devoted to the derivation of the natural orbital and natural 
occupation numbers through the diagonalization of the one-body density matrix. 

In recent years information-theoretic methods play an increasing role for the 
study of quantum mechanical systems. An example is the application of the Max- 
imum Entropy Principle j22j (MEP) to the calculation of the wave function in a 
potential |2H| using as constraints expectation values of simple observables and 
reconstructing a quantum wave function from a limited set of expectation values. 
The idea behind MEP is to choose the least biased result, compatible with the 
constraints of the problem. Thus the MEP provides the least biased description 
consistent with the available relevant information. This is done by employing a 
suitably defined quantum entropy that measures the lack of information associated 
with the distribution of a quantum state over a given known basis. 

Information entropy is important for the study of quantum mechanical systems 
in two cases: first in the clarification of fundamental concepts of quantum mechan- 
ics and second in the synthesis of probability densities in position and momentum 
space |29] . 

In the present work special effort is devoted for the calculation of various quan- 
tum information properties including Shannon entropy, Onicescu energy, Kullback- 
Leibler relative entropy and also Jensen-Shannon divergence. The information 
properties are calculated for the interatomic correlations. In addition the Gross- 
Piatevskii equation is solved giving the density and momentum distribution which 
are employed to calculate the above quantum information properties of the Bose 
gas |26| . The results are compared with those taken in a phenomenological 
way in the framework of the Jastrow correlations. 

The plan of the paper is the following: In Sec. |21the general definitions related 
to the density matrices of a Bose system are considered. Details of the lowest- 
order cluster expansion, analytical expressions and numerical results are reported 
in Sec. |31 In Sec. |1] formulas for the quantum information properties (both for 
the one-and two-body density matrices) are reviewed and analytical results are 
presented. Quantum information properties based on Gross-Piatevskii equation 
are presented in Sec. |21 while the summary of the work is given in Sec. El 
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2 Definition of Density Matrices 



Let ^(ri, r2, . . . , r^) be the wave function describing the trapped Bose gases. In 
the case where this system is composed of non-interacting bosonic atoms at zero 
temperature, all atoms occupy the same single-particle ground state. The many 
body ground state wave function ^'o(ri, r2, . . . , r^) is then a product of A identical 
single particle ground state wave functions. This ground state wave function is 
therefore called the condensate wave function or macroscopic wave function and 
has the form |21j 

^o(ri,r2, • • • ,ryi) = Vo(ri)'(/'o(r2) • • • V'o(rA), (1) 

where ipo{r) is the normalized to one ground-state single-particle wave function 
describing a bosonic atom. It is worth to indicate that Eq. is valid even when 
weak interactions are included. In this case the wave function ^'o(ri, r2, . . . , r^) is 
still, to a very good approximation, a product of A single particle wave functions 
obtained now from the solution of a non-linear Schrodinger equation, the well 
known Gross-Pitaevskii equation. However, in the general case where interactions 
between atoms are included, the ground state wave function ^^(ri, r2, . . . , r^) is 
modified from the simple form of Eq. (^. In that case a percentage of atoms is 
moving from the condensate orbit ipo to higher orbits. 

In the present work we adopt the following normalization of the wave function 
^(ri,r2,...,rA), 

j ^*(ri,r2,...,rA)^'(ri,r2,...,rA)dridr2---(irA = 1, (2) 

where the integration is carried out over the radius vectors ri, r2, . . ., r^. 

A quantity characterizing very important aspects of a Bose gas (as well a 
variety of quantum many-body systems) is the one-body density matrix defined 
as in ^ 

p(ri,r;) = j ^*{ri,r2,...,rA)^{r'i,r2,...,rA)dr2---drA. (3) 

The one-body density matrix is connected to the position- and momentum- 
space properties of the Bose gas and in addition it is the quantity which gives the 
percentage of the condensate of the system. 

The two-body density matrix is a generalization of the one-body density matrix 
and is defined as 

p(ri,r2;ri,r2) = / ^*(ri, r2, rs, • • • , rA)^'(ri, r2, ra, • • • , rA)dr3 • • • drA- (4) 
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The above density matrices are related by the fohowing equation 



P(ri,ri) = j p(ri,r2;ri,r2)(ir2. (5) 

The two-body density matrix is related directly to the interatomic interac- 
tion and its diagonal part provides the two-body density distribution p(ri,r2) 
(expresses the joint probability of finding two atoms at the positions ri and r2, 
respectively), a key quantity of the present work 

/9(ri,r2) = /)(ri,r2;r;,r2) |r;=ri,rJ,=r2 • (6) 

On the other hand the diagonal part of the one-body density matrix is just 
the density distribution of the Bose gas and expresses the probability of finding 
an atom at position ri 

/)(ri) = p(ri,r;)|r,=r;. (7) 
The quantities p{ri) and p(ri,r2) are also related by the following integral 

p(ri) = J p(ri,r2)(ir2. (8) 

Very interesting is also the description of the Bose gas in momentum-space via 
the quantities of the one- and two-body momentum distributions. The two-body 
momentum distribution n(ki,k2) expresses the joint probability of finding two 
atoms with momentum ki and k2 respectively and is given by a particular Fourier 
transform of the corresponding two-body density matrix p(ri, r2; r'j^, r'j) 



n(ki,k2) = j P(ri,r2;ri,r2)exp[iki(ri-ri)]exp[ik2(r2-r2)]dridr'idr2dr2. 

(9) 

The one-body momentum distribution (or simply momentum distribution) 
n(k), expresses the probability of finding an atom with momentum k, and it is 
given by a particular Fourier transform of the one-body density matrix p{ri,r'^) 

nOi) = j p(ri, r'l) exp [ik(ri - r[)] dridr[. (10) 

It can be shown easily that in the case where the Bose gas is described by the 
wave function of Eq. the two-body density matrix is given by 

po(ri,r2;ri,r2) = po(ri, r'i)po(r2, r's), (11) 
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where 

Po(ri,r;)=VS(ri)V'o(r'i). (12) 
From Eqs. ©, and we get 

no(ki,k2) = j^-^ j po{ri,r[)exp\iki{ri - r[)] dridr[ 

^(2^ / ^o(r2,r2)exp [ik2(r2 - 4)] dr2dr2 (13) 
= no(ki)no(k2). 
Finally form Eqs. ((TUl) and ((T^ we get 

"'o(k) = J^^:^^ j ^S(ri)exp [zkri](iri 

(2^ / ^o(r'i)exp [^kr'i] dr'i (14) 

From Eq. (fTl|) it is obvious that '(/'o(^) particular Fourier transform of 

the single particle wave function V'o(r)- 

2.1 Static Structure Factor 

Spectroscopic studies have been used to assemble a complete understanding of the 
structure of atoms and simple molecules j^. The static structure factor S{k) is 
a fundamental quantity, connected with the atomic structure, and is the Fourier 
transform of the radial distribution function g{r). S{k) gives the magnitude of 
the density fluctuation in the system (atomic, molecular, electronic or nuclear) at 
wavelength 27r/A;, where k is the momentum transfer. In recent papers, the Bragg 
spectroscopic method was used to measure S{k) either in the phonon regime 
or/and in the single-particle regime |32j. 

The static structure factor in a finite system is defined as [33] 

5(k) = 1 + ^ / e^'^(^i-'2) [p(ri, ra) - p(ri)p(r2)] dvidv^. (15) 

In the most general case the two-body density distribution p(ri,r2) and the 
one-body density distribution p(r) are connected via the following relation 

p(ri, r2) = Cp(ri)p(r2)/'(ri2) = C p{v^)p{v2)g{n2), (16) 
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where g{r 12) is the radial distribution function and C is the normahzation factor 
which ensures that 

J p{ri,r2)dridr2 = N{N - 1). (17) 

We also consider that 

J p(ri)dri = N, (18) 

where A'^ is the number of the atoms of the Bose condensate. 

In the uncorrelated case (non-interacting gas) the radial distribution function 
is gi^u) = 1 (absence of correlations), and the two-body density distribution 
becomes 

N -I 

p(ri,r2) = ^^p(ri)p(r2). (19) 
Using Eq. ((T^ . Eq. l(T3|) is written as 

5(k) = 1 + 1 y exp [fk(ri - r2)] p{ri)p{r2)[Cg{n2) - l]dridr2. (20) 

Conditions ^ and ^ ensure that S{0) = 0. 

The integration in Eq. H20|) can be performed if the function g{r) is known. 
g{r) must obey the rules g{r = 0) = and lim g{r) 1. The first rule introduces 

r— >cxD 

the repulsive correlations between the atoms and the second the absence of such 
correlations in long distances. In general the form of S{k) is affected appreciably 
from the form of g{r). More specifically, the long range behavior of g{r) affects 
S{k) for small values of k while its short range behavior affects S{k) for large 
values of A; as a direct consequence of the Fourier transform theory. 



2.2 Natural Orbitals and Natural Occupation Num- 
bers 

In the case of the inclusion of the inter-particle interactions between the atoms, 
which give rise to the depletion of the condensate, the one-body density matrix is 
written Plj 

P(ri,ri) = "'oV'o(ri)^o(r'i) +^niV'*(ri)V'j(r'i), (21) 
where rij = 1. The sum niip* {ri)ipi{r'^) is the contribution arising from the 

i i^O 

atoms out of the condensate. The eigenfunctions ipi{r), which are called natural or- 
bitals (NO's), and the eigenvalues n^, called natural occupation numbers (NON's), 
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are obtained by diagonalizing the one-body density matrix through the eigenvalue 
equation 

j p{ri,r[)ipi{r[)dr[ = niipiiri). (22) 

The condition, generally adopted, for the existence of condensation is that there 
should be one eigenvalue rii which is of the order of the number of the particles in 
the trap. 

The NO's ipi{ri) and the NON's rij are obtained by diagonalizing the one- 
body density matrix through the eigenvalue equation (|22)) by expanding first the 
one-body density matrix in a series of Legendre polynomials Pi{x) 

oo 

p{r,r') = p{r,r', cos UJrr') = pi{r,r')Pi{cosLJrT'), (23) 

where pi{r,r') are the coefficients of the expansion 
2/ + 1 

Pl{r,r') = / p{r,r' , cos LUrr') Pl{cOSiOrr') 'i{cosuJrr') ■ (24) 

From the Eqs. (|22j) . (|23|) and (|24jl the eigenvalue equation is written 

nOO 

'^,„A^O/' '^ '2 1 / _ NO, NO 



4n / p,{r, r')^Z''ir')r"dr' = (r), (25) 

Jo 

where ^^i^ir) is the radial part of ipi{r) (V'j(r) = Vni^ ir)Yimi^r)) ■ 

3 Jastrow type Correlated Properties of a Trapped 
Bose Gas 

3.1 Correlated Density Matrices 

A dilute trapped Bose gas can be studied using the lowest-order approximation 
jl3j . In this approximation the two-body density matrix has the form [23. 24^ 

p(ri,r2;r;,r'2)=iVopo(ri,r;)po(r2,r'2)/(|r;-r'2|)/(|ri-r2|), (26) 

where /(|ri — r2|) is the Jastrow correlation function, which depends on the inter- 
particle distance and A'^o is the normalization factor which ensures that 



j P(ri,r2;r'i,r2)|(ri=r;,r2=r^) dridr2 = 1. 
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The diagonal part of p(ri, r2; r'^^, that is the two body density distribution 
takes the form 

p(ri,r2) = iVo/5o(ri)po(r2)/'(r-i2), (27) 
while the one-body density matrix is given by the integral 

p(ri,r;)=iVopo(ri,r;) J /5o(r2, r2)/(|r; - r2|)/(|ri - r2|)dr2. (28) 

The density distribution, which is the diagonal part of p(ri,r'^), can also be 
obtained from the integral 



/o(r) = y p(r,r2)dr2 = iVo/Oo(r) j /3o(r2)/^(|r - r2|)dr2. 



(29) 



In the present work we consider that the atoms are confined in an isotropic HO 
well where the normalized to 1 ground state single particle wave function 'tpo{r) 
has the form of a Gaussian given by the formula 



^o(r) 









262 



where the width b 



-) 

muj I 



1/2 



while the density distribution has the form Pq{v) = |V'o(r)P- The correlation 
function /(ri2) is taken to be of the form 



/(r) = 1 - exp 



W_ 

' 62 



21 



(30) 



where r = |ri — r2|. The correlation function /(r) goes to 1 for large values of r 
and goes to for r — > 0. It is obvious that the effect of the correlations introduced 
by the function f{r), becomes large when the correlation parameter y becomes 
small and vice versa. 

The above defined correlation function was used in [231 124j to find analytical 
expressions of the one-body density matrices in position and momentum spaces 
and static structure factor, while the NO's and NON's are calculated numerically 
employing Eq. (|25|) . 

The analytical expression of the two-body density matrix obtained from Eq. 
(1261) has the form 



/)(ri,r2,ri,r2) 



7^356 



Oi(ri,r2,r'i,r2) - 02i(ri, r2, r'l, rg) 



-022(ri,r2,r'i,r2) + 023(ri, r2, r'^, 



(31) 
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where A'o is the normahzation factor of the form |^ 

2 1 



iVn 



1 



+ 



(l + 22/)3/2 (l+4y)3/2 



(32) 



and 

Oi(ri,r2,r;,r'2) 
02i(ri,r2,r'i,r2) 
022(ri,r2,r'^,r'2) 
023(ri,r2,r;,r2) 



exp 



exp 



exp 



exp 



nb + 4b 



'9 '9 
^ll + Wb 



exp 



exp 



(l + 2y)(ri| + r 



2bJ 



{l + 2y){rl, + r[l + 4, + 4l) 



exp [2yr;^r'2j 
exp [2yrif,r2f,] 
exp [2y(ri5r2fe + r'^^r 



where r^, = r/b. 

The two-body density distribution in accordance with Eq. © is given by 



P(ri,r2) 



exp[-rfj exp[-r|f,] (l - exp[-y(rib - r26)^]^ 



(33) 



The analytical expressions of the one-body density matrix obtained from Eq. 
((2H1) has the form 



P(r,r') 



VT' 



No 
3/253 



[Oi(r, r') - 02i(r, r') - 022(r, r') + 023(r, r')] , (34) 



where the one- and the two-body terms of the expansion in the low order approx- 
imation have the forms 



Oi(r,r') 

02i(r,r') 
022(r,r') 

023(r,r') 



exp 



2 ^ 



exp 



(l+y)3/2 

02i(r',r), 

1 

(l + 2y)3/2 

X exp 



l + 3yrl_rf 
1 + y 2 2 



(35) 

(36) 
(37) 



./2' 



-(l + 2?y)ii^^ 



y ( , f\2 



(38) 
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The analytical expression of the density distribution can be found from Eq. 
(jSH), putting r' = r 



p{r) 



No 
^3/253 

+ 



exp [- 



1 + 2y)3/2 



exp 



(l + y)3/2 



exp 



1 + 27/ 2 



(39) 



The two-body momentum distribution is calculated form the integral of Eq. 
© and has the form 



n(ki,k2) 



66 



X I 1 



Noexp[-klf,] exp[-A;^f,] 



:i + 4y)3/2 



exp 



l + 4y' 



^ib 



k2fe)' 



(40) 



while the momentum distribution can be found analytically using Eq. pOj) and 
has the form 



n{k) 



exp [-k^\ 



^3/2 
+ 



(1 + 3y)3/2 



exp 



1 + 2^ ,2 
l + 3y ^ 



(l + 2y)3/2(l + 4y)3/2 



exp 



1 - 
l + 2y' 



(41) 



where fc;, = kb. 

The above analytical expressions of p{r) and n(A;) have been used to find the 
analytical expressions of the mean square radius and kinetic energy of the trapped 
gas. The expressions we found, for (r^) and (T), are 



and 



(r^) = Nob'' 



(T) = Nohuj 



1 + 2/ 



3 1 + 2y 

2 "(l + 2y)5/2 ' 2(1+ 4y)5/2 



+ 



3 3 l + 3y 3 l + 2y 
4~ 2(l + 2y)5/2 + 4(l + 4y)3/2 



(42) 



(43) 



These expressions, which for a given HO trap are functions of the correlation 
parameter y, could be used to find the value of y from Eq. (|l2|) . if the rms radius 
of the trapped atoms is known and then to define (T) from Eq. (|43j) and vice 
versa. For very large values of y Eqs. H42() and H43|) give the HO expressions of 
(r^) and (T), i.e. (r^) = and (T) = jhuj, respectively. 

The calculation of the density distribution of a trapped Bose gas, confined 
in an isotropic HO potential with length b = 10^ A, has been carried out on 
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the basis of Eq. H39|) j23j . The dependence of the density distribution on the 
parameter y, including also the uncorrelated case (y = oo), has been plotted in 
Fig. n^a). It is seen that, the large values oiy {y > 10) correspond to the Gaussian 
distribution (HO case), while when y becomes small enough [y < 1) the density 
distribution spreads out as in Gross-Pitaevskii's theory j23| . For y > 10 the effect 
of correlations is small, while for very large correlations (y < 0.1) the density 
distribution is modified entirely compared to the Gaussian form originating from 
the HO trap. 




Figure 1: (a) The density distribution p{r) versus = r/b {h = 10* A) for 
various values of the parameter y. The normalization is J p{r)dr = 1. (b) 
The rms radius, \^{rf}, versus the parameter y. 

The rms radius of the Bose gas has also been calculated analytically from Eq. 
for various values of the parameter y. From Fig. ^b), it is seen that ^J (r^) 
(rb = r/6) is a decreasing function of the parameter y and for y > 10 approaches 
the rms radius of the uncorrelated case. 



3.2 Natural Orbitals and Natural Occupations Num- 
bers 

The NO's and the NON's were calculated by diagonalizing the one-body 

density matrix through Eq. H25|) . The NON nis, gives directly the condensation 
fraction no as a result of the repulsive interaction between the atoms of the Bose 
gas at zero temperature. The NON's for the Is, Ip, Id and If states are given 
in Table ^ It seems that, for strong correlations, a fraction of atoms spread 
out into many states. The condensation fraction nis, versus the parameter ^ is 
plotted in Fig. [3 From that figure and from Table ^ it is seen that the effect 
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Figure 2: The condensate fraction riis, at zero temperature, for interacting 
atoms versus 1/y. 
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Table 1: The natural occupation numbers for various values of the correlation 
parameter y [2^1 • 

of the correlations on ni^ is small and all the atoms occupy the Is ground state, 
when y > 10. The effect of the correlations is prominent when y < 10, while the 
decrease of the parameter y (strong correlations) induces a significant depletion of 
the condensated atoms spreading them into many states. 

The NO's of the states Is, Ip and Id for y = 1 are shown in Fig. |3I It is 
seen that the interatomic correlations in the Is-state spread out the ground state 
wave function and consequently the condensation appears in the outer region of 
the trap. From the same figure it is obvious that the NO's of the Ip and Id states 
are much more localized in coordinate space than the equivalent HO orbitals. 

The momentum distribution can be calculated analytically from Eq. (|41|) or by 
Fourier transform of the NO's. The momentum distribution calculated analytically 
for various values of the parameter y has been plotted in Fig. |^a) . It is seen that 
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Figure 3: The NO's (dot lines) versus rf, (rft = r /b) of the states (a) Is, (b) 
Ip, and (c) Id obtained by diagonahzation of the one-body density matrix 
(for y=l). The sohd hues correspond to the HO wave-function with the trap 
length h = 10^ A. 



the large values of y {y > 10) correspond to the Gaussian distribution, while 
when y becomes small enough (y < 1) the momentum distribution has a sharp 
maximum for k = 0. The momentum distribution of the Is NO state as well as 
of the rest of the NO states for y = 0.1 are shown and compared with the total 
momentum distribution in Fig. Efb). It is obvious that although the Is NO state 
gives the main contribution to the momentum distribution, the additional NO 
states contribute to the momentum distribution mainly in the large values of the 
momentum k. 

The dependence of the mean kinetic energy (T) on the parameter y calculated 
analytically, using Eq. H43() . is presented in Fig. E^c). It is seen that (T) has 
a maximum for y ~ 2.5. It is interesting to note that for the same value of the 
parameter y the NON's of the states Id and If have the same value as can be 
seen from Tabled The contribution of the Is NO state and of the rest of the NO 
states to (T) are shown in the same figure. It is seen that, for large values of the 
parameter y (weak correlations) the main contribution to (T) comes from the Is 
NO state, while for strong correlations there is a significant contribution coming 
from the NO's of the additional states. 

A few comments are appropriate. In this section we study the behavior of 
various condensate quantities treated in the Jastrow manner, which introduces 
one parameter. The determination of that parameter could be made by fit of the 
theoretically calculated quantities (density distribution, momentum distribution, 
(r^), and (T)) to the experimental ones as we mentioned in the end of subsection 
13.11 provided that there are experimental data for the corresponding quantities. 
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Figure 4: (a) The normalized to 1 momentum distribution versus kfy [kf, = 
kb) for various values of the correlation parameter y. (b) The momentum 
distribution and the contribution to it from the NO's of the Is NO state and 
of the rest of the NO states (for y = 0.1). (c) The mean kinetic energy per 
atom, {Th) {Th = T /{hjj)), versus y. The solid curve corresponds to the total 
values of (T;,), while dashed and dotted lines are the contributions to the 
total (Tft) of the NO's of the Is NO state and of the rest of the NO states, 
respectively. 

It could be determined also by using the density distribution or the two-body 
density matrix as a trial one and applying the variation principle to the ground 
state energy of the system. The present approach is quite frequent in the study of 
the quantum many body problem when the solution of the Schrodinger equation 
is very difficult. It should be noted also that in the present work there is not a 
direct dependence between the condensation and the number of the atoms. The 
inter-particle correlations are incorporated in the mean field only by the corre- 
lation function which, in some way, depends on the effective size of the atoms. 
That dependence can be found from the information entropy S using the linear 
dependence of S on ln(iVa;,) and the linear dependence of S on ln(i) (see Sec. 

Ell). ^ 
3.3 Static Structure Factor 

In order to calculate the static structure factor in the framework of the atomic 
calculations we choose two trial forms for g{r) The first one is a gaussian 

type which has been extensively and successfully used for the study of similar 
problems in atomic physics (Bose gas, liquid helium) as well in nuclear physics. 
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The relevant g{r) and the entailed ^(A;) (Case 1) are 
g(r) = 1 -exp[-/3r2], 
S{k) = 1 + N{Ci - 1) exp 
NCi 



1.21 



(1 + 2y2)3/2 



exp 



A-2 



l + 2y) 



(44) 



where fc^ = kb, y = pb'^, (5 is the correlation parameter and Ci is the normalization 
factor. 




Figure 5: The static structure factor S{k) of the trapped Bose gas in various 
cases versus the momentum k, (a) in Case 1 for various values of the corre- 
lation parameter j3 as well as for the uncorrelated case (harmonic oscillator) , 
(b) in Case 2 for the least squares best fit value of the parameter a. The 
experimental points are from reference For the various cases see text. 



The second trial function g{r) and the relevant S{k) (Case 2) are of the form 



9{r) 
S{k) 



1 



sin^ ar 



(ar)4 ' 
1 + N{C2 - 1) exp 

5 



a,; 



2 



Aexp 



4 



.(45) 



where a is the correlation parameter, are known coefficients, f3i = (3i{a,b,kb), 
erffz) = —= / e~* dt and is the normalization factor 1251 . 
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Figure 6: The radial distribution function g{r) for Case 1 and 2 (correspond- 
ing to inhomogeneous Bose gas) with the best fit values of the correlation 
parameters. 

The behavior of S(k), in Case 1, for various values of the correlation parameter 
P is shown in Fig. Eta). It is obvious that the effect of correlations, induced by 
the function g{r), becomes large when the parameter /? becomes small and vice 
versa. The case where /3 — > oo, corresponds to the uncorrelated case (HO). For the 
values of k employed in the experiment of Ref. (hereafter EXP) the prediction 
of the HO model is always close to 1 for S{k). When the correlation parameter 
P decreases considerably (strong correlations) the theoretical prediction of ^(A;) is 
in good agreement with the experimental data. The value /3=5.3 fim~^ gives the 
best least squares fit in that case. In general the gaussian form of g{r), in spite 
of its simplicity, reproduces fairly well the experimental data of EXP, both in low 
and high values of the momentum k. Within our theoretical model, the gaussian 
type of g{r) is flexible enough to obtain values for S{k) in agreement with the 
experimental data. 

Fig. Efb) displays the results in Case 2, which are compared with those of the 
data of EXP. The model reproduces well the experimental data in the range 1.5-3 
^m~^ (with best least squares fit value a = 1.34 fj,m~^), but fails in the range 
k > 3 fim^^. The main drawback of this model is the predicted negative values 
of S{k) in the range close to A; = when the correlation parameter a decreases 
considerably (strong correlation case). 

The correlation function g{r) corresponding to Cases 1 and 2 for the correlation 
parameters /?=5.3 fim^^ and a = 1.34 ^m~^ respectively is sketched in Fig. El 
Those values of the parameters (3 and a give the best value in the fit of the 
theoretical expressions of S{k) to the data of EXP. The most striking feature 
in Case 2 is the existence of strong correlations, introduced by g{r)^ in order to 
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reproduce the experimental data of ^(A;). It is worthwhile to point out that g{r), 
in Case 2, exhibits fluctuations in the range r > 2 ^um but this is not visible in 

Fig. m 

The possibility of a linear dependence of S{k) on k for small values of k, as 
predicted from other works j33j . is prohibitive, on the basis of Eq. (|15() at least in 
the case where the trap is an harmonic oscillator one. That can be seen considering 
the ground state wave function to be the harmonic oscillator one and transforming 
ri and r2 in Eq. p5|) into the coordinates of the relative motion (r = ri — r2) and 
the center of mass motion (R = (ri + r2)/2). After some algebra S(k) takes the 
form 

5(k) ~ j exp [ikr] exp [-r^] [Cg{r) - l]dr. (46) 

For finite systems, as a trapped Bose gas, we can expand the exponential 
exp [ikr] , since r is bounded. Thus: 

(ikr)^ (ikr)^ , , 

exp [ikr] = 1 + ?kr + + + • • • . (47) 

Substituting Eq. (|^7)) into Eq. and taking into account that the terms 

with odd powers of k do not contribute to the integral, >S'(A;) takes the form 

S{k) ~ aifc2 + aa/fc^ + • • • . (48) 

Hence, for small values of k, S{k) depends linearly on k'^. The gaussian factor 
exp [— r^] , originating from the harmonic oscillator wave function of the trapped 
Bose gas, ensures the convergence of the integrals Cj corresponding to the even 
powers of the expansion |25j . 



4 Quantum- Information properties of trapped 
Bose gas 

4.1 Shannon Information Entropy 

The Boltzmann-Gibbs-Shannon information entropy |36M37j of a finite probability 
distribution {pi,P2,-iPk) is defined as the quantity 

k 

S = -Y,Pilnpi, (49) 

i=l 

k 

with the constraint: ^^Pi = 1- S is measured in bits if the base of the logarithm 

i=l 

is 2 and nats (natural units of information) if the logarithm is natural. 
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S appears in different areas: information theory, ergodic theory and statistical 
mechanics. It is closely related to the entropy and disorder in thermodynamics. 
The maximum value of S is obtained ii Pi = P2 = ■ ■ ■ = Pk = ^ i-e. Smax = In k. 
The minimum value of S is found when one of the pi's is equal 1 and all the others 
are equal to 0. Then Smin = 0. The above definition holds for discrete probability 
distributions j35j . In quantum mechanics we are often interested in a continuous 
probability distribution p{x). In this case the obvious generalization of Eq. H49|l 
is the information entropy 



where / p{x)dx = 1. Now p{x) is a quantum mechanical probability distribution 
and S may be called the quantum entropy [HHI- "S* indicates the amount of dis- 
order or randomness (uncertainty) in a physical system. Shannon considered this 
uncertainty attached to the system as the amount of information carried by the 
system. If a physical system has a large uncertainty and one obtains information 
on the system by some procedure, as a measurement, then the information is more 
valuable than that received from a system having less uncertainty. Thus, before a 
measurement, the uncertainty of the position of a particle is small for a localized 
probability distribution, while for a diffuse distribution is large. The same holds 
for the missing information due to a limited knowledge of the system via a proba- 
bility distribution. After the measurement the gain in information for a localized 
distribution is smaller than the corresponding gain for a diffuse distribution. 

An important step is the discovery of an entropic uncertainty relation (EUR) 
jH9) ■ which for a three-dimensional system has the form 



where Sr is the information entropy in position-space of the density distribution 
p{r) of a quantum system 



and Sk is the information entropy in momentum-space of the corresponding mo- 
mentum distribution n(k) 




(50) 




(51) 




(52) 




(53) 



The total information entropy is given by 



S — Sf -\- S]^ . 



(54) 
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The density distributions p{r) and n(k) are normalized to one. Inequality 
(jSlj) . for the information entropy sum in conjugate spaces, is a joint measure of 
uncertainty of a quantum mechanical distribution, since a highly localized p{r) 
is associated with a diffuse n(k), leading to low Sr and high Sk and vice- versa. 
Expression (|51|) is an information-theoretical relation stronger than Heisenberg's. 

In previous work we proposed a universal property of S for the density distri- 
butions of nuclei, electrons in atoms and valence electrons in atomic clusters |4U| . 
This property has the form 

S = a + blnN, (55) 

where is the number of particles of the system and the parameters a, b depend 
on the system under consideration. It is noted that recently we have obtained 
the same form for systems of correlated bosons in a trap [201 • This concept was 
also found to be useful in a different context. Using the formalism in phase- 
space of Ghosh, Berkowitz and Parr |41j . we found that the larger the information 
entropy, the better the quality of the nuclear density distribution [l^l- Recently 
the Shannon information entropy has been applied successfully to the study of the 
free expansion of impenetrable bosons on the one-dimensional optical lattices jl3] • 



4.2 Onicescu's Information Entropy 



Onicescu tried to define a finer measure of dispersion distributions than that of 
Shannon's information entropy 01]. Thus, he introduced the concept of informa- 
tion energy E. For a discrete probability distribution {pi,P2, ■ ■ ■ ,Pk) the informa- 
tion energy E is defined by 



i 

which is extended for a continuous density distribution p{x) as 



E 



p^{x) dx. 



(56) 



(57) 



The meaning of (|57|) can be seen by the following simple argument: For a Gaussian 
distribution of mean value p, standard deviation a and normalized density 



p{x) 



1 



27rcr 



exp 



2ct2 



relation gives 



E 



27rcr2 



poo 


{x-pf 


/ exp 




J —oo 


Cj2 



dx 



2cr/7r' 



(58) 



(59) 
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E is maximum if one of the pj's equals 1 and all the others are equal to zero i.e. 
Emax = 1, while E is minimum when pi = p2 = . . . = pk = \, hence Emin = \ 
(total disorder). E has been called information energy, although it does not have 
the dimension of energy j45j . This is due to the fact that E becomes minimum for 
equal probabilities (total disorder), by analogy with thermodynamics. 

It is seen from (|59|) that the greater the information energy, the more concen- 
trated is the probability distribution, while the information content decreases. E 
and information content are reciprocal, hence one can define the quantity ^24^ 

= ^, (60) 

as a measure of the information content of a quantum system corresponding to 
Onicescu's information energy. 

Relation (|57jl is extended for a 3-dimensional spherically symmetric density 
and momentum distribution as follow p(r) 

Er= p^{r)dT 



Ek = j n2(k)dk. (61) 

Er has dimension of inverse volume, while E^ of volume. Thus the product 
ErEk is dimensionless and can serve as a measure of concentration (or information 
content) of a quantum system. It is also seen from (|59)) . (|60|) that E increases 
as a decreases (or concentration increases) and the information (or uncertainty) 
decreases. Thus O and E are reciprocal. In order to be able to compare O with 
Shannon's entropy S, we redifine O as 

O = (62) 

as a measure of the information content of a quantum system in both position and 
momentum spaces, inspired by Onicescu's definition. 

4.3 Landsberg's Order Parameter 

Landsberg defined the order parameter $7 (or disorder A) as 

where S is the information entropy (actual) of the system and 5 (max) the maxi- 
mum entropy accessible to the system. Thus the concepts of entropy and disorder 
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are decoupled and it is possible for the entropy and order to increase simultane- 
ously. It is noted that $7 = 1 corresponds to perfect order and predictability, while 
= means complete disorder and randomness. 



4.4 Two-body information entropies 

The two-body Shannon information entropy both in position- and momentum- 
space and in total are defined respectively 071 UHl 

S2r = - J /9(ri,r2)lnp(ri,r2)(iri(ir2 (64) 
S2k = - J n(ki,k2)lnn(ki,k2)(ikidk2, (65) 

S2 = S2r + S2k- (66) 

The one-body Onicescu's information entropy is already defined in (|61j) and 
H62|) . where the generalization to the two-body information entropy is straightfor- 
ward and is given by 

O2 = (67) 

^2r-t^2k 



where 



E2r = j p^(ri,r2)dridr2 

E2k = j n2(ki,k2)(iki(ik2. (68) 



4.5 Kullback-Leibler relative entropy and Jensen-Shannon 
divergence 

A well known measure of distance of two discrete probability distributions Pi'\pf^ 
is the Kullback-Leibler relative entropy [IH] 

K{pf\pf^) = Y,^p\nP^ (69) 

i Pi 

which for continuous probability distributions defined as 

K = j p(i)(x)ln^^dx, (70) 
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which can be easily extended for 3-dimensional systems. 

Our aim is to calculate the relative entropy (distance) between p^^^ (correlated) 
and p(^) (uncorrelated) densities both at the one- and the two-body levels in order 
to assess the influence of short range correlations (SRC) through the correlation 
parameter y, on the distance K j24j . It is noted that this is done for both systems 
under consideration: nuclei and trapped Bose gases. An alternative definition of 
distance of two probability distributions was introduced by Rao and Lin [201 1^ , 
i.e. a symmetrized version of K, the Jensen-Shannon divergence J j52j 

J(J>'".P'«) = H jEi^J _ }.„ _ l„ (,(2)) , (71) 

where H{p) = — ^^pjlnpj stands for Shannon's entropy. We expect for strong 

i 

SRC the amount of distinguishability of the correlated from the uncorrelated dis- 
tributions is larger than the corresponding one with small SRC. We may also see 
the effect of SRC on the number of trials L needed to distinguish p^^^ and p*-^-* (in 
the sense described in [SH). 

In addition to the above considerations, we connect Sr and Sk with fundamen- 
tal quantities i.e. the root mean square radius and kinetic energy respectively. We 
also argue on the effect of SRC on EUR and we propose a universal relation for 
S, by extending our formalism from the one- and two-body level to the A'^-body 
level, which holds exactly for uncorrelated densities in trapped Bose gas and it is 
conjectured to hold approximately for correlated densities in Bose gases (see Sec. 

EZI). 

The Kullback-Leibler relative information entropy K for continuous distribu- 
tions p[^^ and p[^^ is defined by relation (|7n|). It measures the difference of p[^^ form 

(2) 

the reference (or apriori) distribution . It satisfies: K > for any distributions 
p^^^ and /9^^\ It is a measure which quantifies the distinguishability (or distance) 
of from p^-^^ , employing a well-known concept in standard information theory. 

In other words it describes how close p[^^ is to p[^^ by carrying out observations or 
coin tossing, namely L trials (in the sense described in [HII)- We expect for strong 
SRC the amount of distinguishability of the correlated p[^^ and the uncorrelated 

(2) 

distributions p- is larger than the corresponding one with small SRC. 

However, the distance K does not satisfy the triangle inequality and in addition 
is i) not symmetric ii) unbounded and iii) not always well defined To avoid 
these difficulties Rao and Lin |5n| IBlj introduced a symmetrized version of K 
(recently discussed in |52j^. the Jensen-Shannon divergence J defined by relation 
(|71|) . J is minimum for p^^^ = p^"^^ and maximum when p^^^ and p^'^^ are two 
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distinct distributions, when J = ln2. In our case J can be easily generalized 
for continuous density distributions. For J minimum the two states represented 
by p^^^ and p^"^^ are completely indistinguishable, while for J maximum they are 
completely distinguishable. It is expected that for strong SRC the amount of 
distinguishability can be further examined by using Wooter's criterion |S2]. Two 
probability distributions p^^^ and p^"^^ are distinguishable after L trials (L — > oo) 
if and only if /j^^))) 5 > 

The relative entropy is a measure of distinguishability or distance of two states. 
It is defined generahzing (jTO]). by 

^ = /^'(r)ln|5|^^r. (72) 

In our case ip{r) is the correlated case and (j){r) the uncorrelated one. Thus 

K^r= I P(r) In^T^dr, (73) 

where p{t) is the correlated one-body density and p'{r) is the uncorrelated one- 
body density. 

A corresponding formula holds in momentum-space 

K,, = y"n(k)ln^dk, (74) 

where n(k) is the correlated one-body density and n'(k) is the uncorrelated one. 
For the two-body case we have 

K2r = I p(ri,r2) In^T^il^dndrs, (75) 

where p{ri,T2) is the correlated two-body density in position-space and p'{ti,T2) 
is the uncorrelated one. 

The generalization to momentum-space is straightforward 

K2k = I n(ki, k2) In (76) 
J n'(ki,k2) 

where n(ki, k2) is the correlated two-body density in momentum-space and n'(ki, k2) 
is the uncorrelated one. 

For the Jensen-Shannon divergence J we may write formulas for Ji (one-body) 
and J2 (two-body), employing definition (|71l) and putting the corresponding corre- 
lated p*-^-* and uncorrelated p^"^^ distributions in position- and momentum-spaces. 
We calculate K and J in position- and momentum-spaces, for nuclei and bosons. 
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4.6 Numerical Results and Discussion 



For the sake of symmetry and simplicity we put the width of the HO potential 
6 = 1. Actually for 6 = 1 in the case of uncorrelated case it is easy to see that 
Sir = Sik and also S2r = S2k (the same holds for Onicescu entropy), while when 
6/1 there is a shift of the values of Sir and Sik by an additive factor In 6^. 
However, the value of 6 does not affect directly the total information entropy S 
(and also O). S and O are just functions of the correlation parameter y |24j . 

In Fig. [7|we present the Shannon information entropy Si using relation H54|) 
and 5*2 using relation (|66|) in trapped Bose gas as functions of the correlation 
parameter lii(^). It is seen that 5*1 and 5*2 increase almost linearly with the 
strength of SRC i.e. In (i) in both systems. The relations 52 = and O2 = 0\ 
hold exactly for the uncorrelated densities while the above relations are almost 
exact for the uncorrelated densities. For the sake of comparison we also present 
the decomposition of S in coordinate and momentum spaces i.e. Sir-, Sik-, S2r, S2k 
employing (|52|) . (|53|) . (|65j) . The most striking feature concluded from the 

above Figures is the similar behavior between Sir and S2r and also Sik and S2k 
respectively. 




Figure 7: (a) Shannon information entropy (one- and two-body) (b) Shan- 
non information entropy (one- and two-body) both in coordinate- and 
momentum-space, in a trapped Bose gas. 

In Fig. IS] we plot the Onicescu information entropy both one-body (Oi) and 
two-body (O2) (relations (jHTjl ). We conclude by noting once again the strong 
similarities of the behavior between one- and two-body Onicescu entropy. 

It is interesting to observe the relation of the rms radii \/Jr^ with Sr as well 
as the corresponding relation of the mean kinetic energy (T) with Sk, as functions 
of the strength of SRC, In(^). This is done in Fig. IHl for (r^) and (T) after 
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Figure 8: The Onicescu information entropy (both one- and two-body) in a 
trapped Bose gas. 

applying the suitable rescaling. The corresponding curves are similar for nuclei 
and trapped Bose gas. 



3.50 



3.25 



"""-4 -2 """-4 -2 

ln{1/y) ln(1/y) 

Figure 9: (a) Mean-square radius and the Shannon information entropy 
Sir (b) Mean kinetic energy (T) (in hw units) and the Shannon information 
entropy Sik, as functions of the correlation parameter In (^) in a trapped 
Bose gas. 

A well-known concept in information theory is the distance between the prob- 
ability distributions p[^^ and p^'^\ in our case the correlated and the uncorrelated 
distributions respectively. A measure of distance is the Kullback-Leibler relative 
entropy K defined previously. The correlated and uncorrelated cases are compared 
for the one-body case (i^i) and for the two-body case {K2) in Fig. IIOI decompos- 
ing in position- and momentum-spaces according to (|73|) - (|76j) . It is seen that Kir, 
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increase as the strength of SRC increases, while i^ifc, have a maximum 
at a certain value of In {-) depending on the system under consideration. 




Figure 10: (a) One- body KuUback-Leibler relative entropy both in 
coordinate- and momentum-space (b) Two-body KuUback-Leibler relative 
entropy both in coordinate- and momentum-space, in a trapped Bose gas. 

Calculations are also carried out for the Jensen-Shannon divergence for one- 
body density distribution ( Ji entropy) as function of In (^), decomposed in position- 
and momentum- spaces (Fig. We observe again that Ji increases with the 

strength of SRC in position-space, while in momentum-space there is a maximum 
for a certain value of In (-) . It is verified that < J < In 2 as expected theoretically 

It is noted that the dependence of the various kinds of information entropy 
on the correlation parameter lii(i) is studied up to the value lii(^) = {y = 1), 
which is already unrealistic corresponding to strong SRC. In addition, lowest order 
approximation does not work well beyond that value. In this case three-body terms 
should be included but this prospect is out of the scope of the present work. 

For very strong SRC the momentum distribution n{k) exhibits a similar be- 
havior with the mean field {y oo). This is illustrated in Fig. 13 where we present 
n{k) for various values of In(^). It is seen that for small and large SRC the tail of 
n{k) disappears. That is why for small and large SRC the relative entropy (i^ifc 
and Jik) is small, while in between shows a maximum (Fig. I1U|1 . A similar trend 
of n(ki, for large SRC explains also the maximum of the relative entropy K2k- 

4.7 Conclusions and comments 

Our main conclusions are the following [2^ 
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Figure 11: (a) One-body Jensen-Shannon divergence entropy both in 
coordinate- and momentum-space, in a trapped Bose gas (b) Momentum dis- 
tribution n{k) of ^He for various values of the correlation parameter lii(^). 
The case MF (mean field) corresponds to the uncorrelated case {y oo). 

(i) Increasing the SRC (i.e. the parameter ln(i)) the information entropies 
S, O, K and J increase. A comparison leads to the conclusion that the 
correlated Bose gas has larger values of entropies than the uncorrelated one. 

(ii) There is a relation of with (r^) and Sj. with (T) in the sense that 
they have the same behavior as a function of the correlation parameter 
ln(i). These results can lead us to relate the theoretical quantities Sr and 
Sk with experimental ones like charge density distribution, and momentum 
distribution, radii, etc. A recent paper addresses this problem |5,Sj . 

(iii) The relations 5*2 = and O2 = Of hold exactly for the uncorrelated 
densities in trapped Bose gas while the above relations are almost exact in 
the case of correlated densities. In previous work we proposed the universal 
relation Si = Sr + Sk = a + b In N where is the number of particles of the 
system either fermionic (nucleus, atom, atomic cluster) or bosonic (corre- 
lated atoms in a trap). Thereby, in the general case (including correlations 
among the atoms) 

^2 ~ 2^1 = 2(a + 6 InA). 

For the generalized A-body uncorrelated distributions /9(ri,r2, . . . ,r]\f) and 
n(ki, k2, . . . , k^v) the relation 

Sn = NSi = A (a + 6 In A). 

holds exactly. 
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It is conjectured that it holds approximately for correlated systems (which 
has still to be proved for > 3). 

(iv) The entropic uncertainty relation (EUR) is 

S = Sr + Sk> 6.434. 

It is well-known that the lower bound is attained for a Gaussian distribution 
(i.e. the case of uncorrelated Bose gas). In all cases studied in the present 
work EUR is verified. 

A final comment seems appropriate. In general, the calculation of p(ri,r2) 
and n(ki, k2) is a problem very hard to be solved in the framework of short 
range correlations. In the present work we tried to treat the problem in 
an approximate but self-consistent way in the sense that the calculations of 
p(ri,r2) and n(ki,k2) are based in the same p(ri, r2; r'^, r2), which is the 
generating function of the above quantities. As a consequence the informa- 
tion entropy S2 = S2r + 5'2fc is derived also in a self-consistent way and there 
is a direct link between S2r and S2k , as well as the other kinds of information 
entropies which are studied in the present work. 



5 Quantum- Information properties based on 
Gross-Pitaevskii equation 

The ground-state properties of the condensate, for weakly interacting atoms, are 
explained quite successfully by the non-linear equation, known as Gross-Pitaevskii 
(GP) equation, of the form 



+ -mu?r^ + N -HMr 

2m 2 m '^^ ^' 



V^(r) =/x^(r), (77) 



where A'^ is the number of the atoms, m is the atomic mass, oq is the scattering 
length of the interaction and /i is the chemical potential This equation has the 
form of a non-linear stationary Schrodinger equation, and it has been solved for 
several types of traps using various numerical methods 1541 The presence 
of the third term, which is linear in N is responsible for the dependence of the gas 
parameter x o^i the density of the system. 

Eq. (f77j) was solved numerically in Ref. for trapped boson-alkali atoms 
in two cases. For a system of ^''Rb atoms with parameter 69 = 12180 A (angular 
frequency ujqI2t: = 77.78 Hz) and scattering length oq = 52.9 A jj^, and for 
a system of ^^^Cs atoms with parameters 69 = 27560 A {ujq /27r = 10 Hz) and 
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ao = 32 A In these cases the effective atomic size is small compared both to 
the trap size and to the interatomic distance ensuring the diluteness of the gas. 

There is in atomic physics a connection of Sr and Sk with the total kinetic en- 
ergy T and mean square radius of the system through rigorous inequalities derived 
using the EUR ^ EHl 



5r(min) < Sr < S'r(max) 
S'fc(min) < 5fc < 5/c(max) (79) 



(78) 

, (79) 
5(min) <S< 5(max). (80) 



The lower and upper limits are written here more conveniently in the following 
form, for density distributions normalized to one: 

5,(min) = ^(l + ln^)-^ln(^r 

5.(max) = ^(l + ln^) + ^lnQ(r2)), (81) 



3.. . , . 3 , /2 , 2, 



5fc(min) = -(l+lnvr)--ln(^-(r^) 
Skimax) = ^(l + ln7r) + ^ln(^^r) , (82) 



S'(min) = 3(l+ln7r) 

5(max) = 3(l+ln7r) + ^ln(^^(r2)ry (83) 

In Ref. [Sni it was verified numerically that the above inequalities hold for 
nuclear density distributions and valence electron distributions in atomic clusters. 
We also found a link of S with the total kinetic energy of the system T, and a 
relationship of Shannon's information entropy in position-space Sr with an exper- 
imental quantity i.e. the rms radius of nuclei and clusters. 

It has been verified numerically [HOI that inequalities (|78|) . (|79|) and (|8fl|) hold 
for correlated bosonic systems as well, i.e. trapped boson-alkali atoms ^''Rb and 
^^^Cs. That is shown in Table [21 for ^''Rb bosonic system. An analogous table 
can be displayed for ^^^Cs. It is noted that for large A^, Sk may become negative, 
but the important quantity is the net information content S = Sr + Sk oi the 
system which is positive. We employed density distributions p{r) and n(k) for 
bosons derived by solving numerically the GP equation ()77() . It is also seen that 
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N 






Sf 


Sk 


Sk 


Sk 


S 


S 


S 


5x10^ 


3.797 


3.834 


3.845 


2.590 


2.630 


2.637 


6.434 


6, 


.465 


6.482 


10 


4.027 


4.100 


A "1 c~\r\ 

4.120 


O O "1 /I 

2.314 


o on A 

2.394 


c\ A r\C) 

2.408 


A ^ A 

6.434 


6, 


.494 


6.528 


3x10^ 


4.437 


4.599 


4.640 


1.794 


1.963 


1.997 


6.434 


6, 


.562 


6.637 


5x10^ 


4.641 


4.855 


4.907 


1.527 


1.746 


1.794 


6.434 


6, 


.601 


6.701 


7x10^ 


4.778 


5.029 


5.090 


1.345 


1.598 


1.657 


6.434 


6, 


.627 


6.746 


10^ 


4.925 


5.219 


5.287 


1.148 


1.437 


1.509 


6.434 


6, 


.655 


6.796 


5x10^ 


5.615 


6.113 


6.211 


0.223 


0.667 


0.819 


6.434 


6, 


.780 


7.030 


10^ 


5.922 


6.511 


6.619 


-0.185 


0.317 


0.512 


6.434 


6, 


.828 


7.132 


5x10^ 


6.654 


7.452 


7.577 


-1.142 


-0.533 


-0.220 


6.434 


6, 


.919 


7.357 


10^ 


6.993 


7.864 


7.992 


-1.557 


-0.920 


-0.560 


6.434 


6, 


.943 


7.432 



Table 2: Values of 5*,., 5*^ and 5* versus the number of particles N for ^'''Rb 
bosonic system (see inequalities (f75|). ^7^. and (|HU|)) ^30] . 



the right-hand-side of inequahty H78|) is nearly an equahty. Thus there exists a 
relation between Sr and T for bosons as well as for nuclei [SH] . 

In Ref. [HU' we addressed the problem of finding Sr and Sk (i.e. the extent 
of p(r) and n(k)) for bosonic many-body systems in order to compare with cor- 
responding results for fermionic systems. First we review the results of Ref. j4Uj 
for systems of fermions where we proposed a universal property for S for the den- 
sity distributions of nucleons in nuclei, electrons in atoms and valence electrons in 
atomic clusters. This property has the form 

S = a + blnN, (84) 

where the parameters a and b depend on the system under consideration. The 
values of the parameters are the following 

a = 5.325, b = 0.858 (nuclei) 

a = 5.891, b = 0.849 (atomic clusters) (85) 
a = 6.257, b = 1.007 (atoms). 

Next we verified (|8l|) employing densities p{r) and n{k) for trapped bosons 
solving the GP equation (|77|) . p{r), derived in this way, and n{k) derived by 
Fourier transform of ip{r), were inserted into equations ()52() and (|53() to find the 
values of Sr, Sk and S = Sr + Sk s^s functions of the number of bosons A^. The 
results are shown in Fig. I12f a). The circles and the triangles correspond to the 
calculated values for the bosonic systems ^^Rb and ^^^Qg^ respectively, while the 
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6 8 10 12 14 2 4 6 8 

In N ln(NaJ 



Figure 12: The information entropy S (a) versus In and (b) versus In(A^ab) 
for the bosonic systems according to the fitted expressions (jMjl and (jHHjl 
respectively. The open circles (*^Rb) and triangles (^^^Cs) come from the 
numerical solution of the GP equation. 



lines to the fitted form of Eq. (|84j) where 

a = 6.0291, b = 0.0678, 5 x 10^ < iV < 10^ (^"^Rb) 

a = 5.9614, b = 0.0661, 10^ < < 5 x 10^ (^^^Cs). (86) 

From the values of the parameters a and 6, for the two systems we examined, 
and from the strength of the interatomic interaction, which is proportional to 
Nao/m we can conclude that for the same trap (same value of bo) and for different 
values of the scattering length or /and of the atomic mass, the linear dependence of 
S" on In does not change. The only change will be in the values of the parameter 
a of Eq. (|84j) . Thus the change of the scattering length or / and of the atomic mass 
will produce a parallel displacement of the lines of Fig. I12f a). From the values 
of the parameter b of Eq. (|84j) which give the slope of the lines of Fig. I12r a) 
corresponding to different sizes of the trap bo and because the GP equation can 
be written in the form 



1^2 1 
•-V + - 
2 2 



'^l + 7;r! + ^7rNab\Airb)\' 



ipbin) = P'b'^bin), (87) 



where = r/bo, ah = ao/b^, fih = n/fkoo and ipbi^^b) = &o V'(r), we should 
conclude that there is a parallel displacement of the lines for the various values of 
bo. That leads us to fit the numerical values of S, for the two systems we examined, 
using the formula 

S = a + bln{Nab). (88) 
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The new values of a and b are now a = 6.3976 and b = 0.0678 for ^'^Rh and 
a = 6.4081 and b = 0.0661 for ^^^Cs. As the two hnes are almost the same we use 
the same parameters a and b for the two bosonic systems 

a = 6.4028, b = 0.0669 (89) 

which are the mean values of the corresponding parameters of the two systems. 
The results are shown in Fig. I12f b'). It is seen that the numerical values of S for 
the two bosonic systems are very close to those calculated from Eq. (|88j) with the 
parameters a and 6 given by Eq. 

The results of S and S'(max), displayed in Table 121 allowed us to calculate the 
order parameter $7 (relation ()63() ) as function of the number of particles in a 
system of trapped correlated boson-alkali atoms. The dependence of on iV is 
shown in Fig. for ^^Rb (an analogous figure can be displayed for ^^^Cs). It 
is seen that Q is an increasing function of N. A similar trend has been observed 
in Fig. 1 of Ref. where Q{N) was calculated for nucleons in nuclei and 

valence electrons in atomic clusters. As stated in j60j . our result is in a way 
counter-intuitive and indicates that as particles are added in a correlated quantum- 
mechanical system, the system becomes more ordered. The authors in jjl, studied 
disorder and complexity in an ideal Fermi gas of electrons. They observed that for 
a small number of electrons the order parameter 0, is small, while increases as 
one pumps electrons into the system and the energy levels fill up. 




Figure 13: The order parameter f2 as a function of the number of particles 
N for ^^Rb bosonic system. 
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5.1 Shannon's Information Entropy Using Density Dis- 
tribution of Correlated Bosons, Results and Dis- 
cussion 

From the analytical expressions of p{r) and n{k) (Eqs. (|39|) and (|41|) respectively) 
for a correlated bosonic system, the information entropy S can be found using 
Eqs. (|52|) and H53I) . The values of 5 as function of the correlation parameter ^ are 
shown in Fig. I14f a). The open squares correspond to the calculated values of S 
while the line to the fitting expression 

S = a + bln{l/y), (90) 

where 

a = 6.6687, b = 0.0913, 0.05 < y < 10. (91) 




Figure 14: (a) The information entropy S versus the correlation parameter 
y for a bosonic system. The line corresponds to the fitted expression (j9Up . 
while the open squares to values found using the JCM. (b) The information 
entropy S versus In(iVab) for the bosonic systems according to the fitted 
expression (jHHj) . The open circles (^''Rb) and triangles (^^^Cs) come from the 
numerical solution of the GP equation, while the open squares correspond to 
the JCM. 

We have fitted our numerical results for 0.05 < y < 10. As mentioned before 
large values oiy {y > 10) correspond to a Gaussian distribution, while for y < 0.05 
(very strong correlations) higher order terms must be included in the expansion of 
the density. 
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It should be noted that in the present approach there is not a direct depen- 
dence between the condensation and the number of the atoms. The inter-particle 
correlations are incorporated in the mean field only by the correlation function 
which, in some way, depends on the effective size of the atoms. We could find that 
dependence making the assumption that the correlated parameter y depends on 
N and a;, through the relation 

^ = {\iNa^,f' (92) 

and try to find Ai and A2 equating the rhs of Eqs. (|88|) and (|90j) . In this way Ai 
and A2 can be found as functions of the parameters a and h of Eqs. ()88() and H90|) 
having the forms 

=e(ai-a2)/6i^ A2 = 61/62 (93) 

where oi and 61 are the values of the parameters a and 6 of Eq. H89|) and 02 and 62 
the parameters of Eq. (|91j) . The numerical values of Ai and A2 are: Ai = 0.0188 
and A2 = 0.7330. 

The values of S for the bosonic systems ^^Rb and ^^^Cs found in Jastrow 
correlation method (JCM), versus In(A^afe) (calculated from Eqs. (|^^ and 
are shown in Fig. I14r b) with open squares. It is seen that the two bosonic systems 
studied with the GP theory and the bosonic system studied with the JCM give 
very similar results for S. It seems that the information entropy S for the bosonic 
systems depends only on In(A^afe). 



6 Summary 

The effect of the interparticle correlations between Bose atoms at zero temperature 
is examined using a phenomenological way to incorporate the atomic correlations. 
This is made by introducing the Jastrow correlation function in the two-body 
density matrix. Analytical expressions are found for the one- and the two-body 
density and momentum distribution, mean-square radius, kinetic energy and static 
structure factor. The introduction of correlations changes the shape of the density 
and momentum distributions compared with the Gaussian one, corresponding to 
the harmonic oscillator model. There is a decrease of the density distribution in 
the central region of the atomic system while the momentum distribution increases 
in the region of small k and thus there is a decrease of the mean kinetic energy of 
the system. In addition the natural orbitals and the natural occupation numbers 
have been calculated and consequently the condensate fraction has been obtained 
for different values of the parameter y. A theoretical calculation of the static 
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structure factor is reported also by applying two trial forms for the radial distribu- 
tion function. Our results are compared with recent experimental data concerning 
trapped Bose gas. By applying suitable parametrization the experimental data 
are reproduced quite well. 

Various kinds of quantum information properties of the trapped Bose gas are 
calculated i.e. the Shannon and Onicescu information measures for the correlated 
and uncorrelated cases which are compared as functions of the strength of the 
short range correlations. It can be seen that increasing the short range correlations 
the information entropies S and O increase. There is a relation between \/Jr^ 
and Sr and between (T) and S^- It is also conjectured that the relation Sn = 
N(a + b InN) holds approximately for the correlated system. The Gross-Pitavskii 
equation is solved in order to calculate the information properties of the Bose gas 
from another point of view. It is concluded that the Shannon information entropy 
obeys the functional form S = a + blnN. Finally it is shown that Landsberg's 
order parameter O is an increasing function of the number of Bose atoms A'^. 

Acknowledgments 

The work of K. Ch. Chatzisawas and C. P. Panos was supported by Herakleitos 

Research Scholarships (21866) of EIIEAEK and the European Union while the 
work of S. E. Massen, Ch. C. Moustakidis and C. P. Panos was supported by the 
Pythagoras II Research project (80861) of EIIEAEK and the European Union. 



36 



References 

[1] S. N. Bose, Z. Phys. 26 (1924) 178. 

[2] A. Einstein, Sitzber. Kgl. Preuss. Akad. Wiss. (1924) 261; Sitzher. Kgl. Preuss. 
Akad. Wiss. (1925) 3. 

[3] Bose Einstein Condensation, edited by A. Griffin, D. W. Snoke, and 
S. Stringari (Cambridge University Press, Cambridge, 1995). 

[4] F. Dalfavo, S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71 
(1999) 463. 

[5] A. J. Leggett, Rev. Mod. Phys. 73 (2001) 307. 

[6] A. S. Parkins, and D. F. Walls, Phys. Rep. 303 (1998) 1. 

[7] A. L. Fetter, | cond-mat/ 981 1366] 

[8] Ph. W. Courteille, V. S. Bagnato, and V. I. Yukalov, Laser Phys.ll (2001) 
659. 

[9] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cor- 
neU, Science 269 (1995) 198. 

[10] G. Baym and C. Pethick, Phys. Rev. Lett. 76 (1996) 6. 

[11] F. Dalfavo and S. Stringari, Phys. Rev. A 53 (1996) 2477. 

[12] B. D. Esry, Phys. Rev. A. 64 (2001) 063610. 

[13] A. Fabrocini and A. Polls, Phys. Rev. A 60 (1999) 2319. 

[14] A. Fabrocini and A. Polls, Phys. Rev. A 64 (2001) 063610. 

[15] J. L. DuBois and H. R. Clyde, Phys. Rev. A 63 (2001) 023602. 

[16] A. Minguzzi, P. Vignolo, and M. P. Tosi, Phys. Rev. A 62 (2000) 023604. 

[17] M. Naraschewski and R. J. Glauber, Phys. Rev. A 59 (1999) 4595. 

[18] L. Pitaevskii and S. Stringari, Phys. Rev. Lett. 83 (1999) 4237. 

[19] L. P. Pitaevskii, Zh. Eksp. Teor. Phys. 40, 646 (1961); Sov. Phys. Jetp 13 
(1961) 451. 

[20] E. P. Gross, Nuovo Cimento 20 (1961) 454. 

37 



[21 

[22; 

[23 
[24 

[25; 
[26 

[27; 

[28 
[29 
[30 

[31 



[32; 

[33 

[34; 
[35; 

[36 

[37; 

[38 



W. Ketterle, D.S. Durfee, and D.M. Stamper-Kurn, |cond-mat /9904034| 
P. O. Lowdin, Phys. Rev. 97 (1955) 1474. 

Ch. C. Moustakidis, and S. E. Massen, Phys. Rev. A 65 (2002) 063616. 

Ch. C. Moustakidis, K. Ch. Chatzisavvas, and C. P. Panos, Int. J. Mod Phys. 
E 14 (2005) 1087. 

Ch. C. Moustakidis, Phys. Lett. A 328 (2004) 319. 

S. E. Massen, Ch. C. Moustakidis, and C. P. Panos, Focus on Boson Research, 
ed. A. V. Ling (Nova, New York, 2005) p. 115. 

J. N. Kapur Maximum- Entropy Models in Science and Engineering., (John 
Wiley, 1989). 

N. Canosa, A. Plastino, R. Rossignoh Nucl. Phys. A 550 (1992) 453. 
P. Garbaczewski, Entropy 7 No4 (2005) 253. 

S. E. Massen, Ch. C. Moustakidis, and C. P. Panos, Phys. Let. A 299 (2002) 
131. 

D. M. Stamper-Kurn, A. P. Chikkatur, A. GorUtz, S. Inouye, S. Gupta, 
D. E. Pritchard and W. Ketterle, Phys. Rev. Lett. 83 (1999) 2876; J. Stenger, 
S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, D. E. Pritchard and W. Ket- 
terle, Phys. Rev. Lett. 82 (1999) 4569. 

J. Steinhauer, R. Ozeri, N. Katz, and N. Davidson, Phys. Rev. Lett. 88 (2002) 
120407. 

F. Zambehi, L. Pitaevskii, D. M. Stamper-Kurn and S. Stringari, Phys. Rev. 
A 61 (2000) 063608. 



S. Stringari cond-mat/0101299 



K. Ch. Chatzisavvas, and C. P. Panos, Int. J. Mod Phys. E 14 (2005) 653. 
C. E. Shannon, Bell Syst. Tech. 27 (1948) 379. ibid. 27 (1948) 623. 
J. J. Halhwen, Phys. Rev. D 48 (1993) 2739. 

M. Ohya, and D. Petz, Quantum Entropy and Its Use ( Springer- Verlag, 
Berhn; New York, 1993). 



[39] I. Bialynicki-Birula, and J. Mycielski, Commun. Math. Phys 44 (1975) 129. 

38 



[40] S. E. Massen, and C. P. Panos, Phys. Lett. A246 (1998) 530. 

[41] S. K. Ghosh, M. Berkowitz, and R. G. Parr, Proc. Natl. Acad. Sc. USA 81 
(1984) 8028. 

[42] G. A. Lalazissis, S. E. Massen, C. P. Panos, and S. S. Dimitrova, Int. J. Mod. 
Phys. E7 (1998) 485. 

[43] M. Rigol, and A. Muramatsu, Mod. Phys. Lett. B 19 (2005) 861. 

[44] O. Oniccscu, R. Acad. Sci. Paris A263 (1996) 25. 

[45] C. Lepadatu, and E. Nitulescu, Acta Chim. Slov. 50 (2003) 539. 

[46] P. T. Landsberg, Phys. Lett. A 102 (1984) 171. 

[47] C. Amovilli, N. H. March, Phys. Rev. A69 (2004) 054302. 

[48] T. M. Cover, and J. A. Thomas, Elements of Information Theory, (Wiley- 
Interscience, New York 1991). 

[49] S. Kullback, Statistics and Information theory, (Wiley, New York, 1959). 

[50] C. Rao, Differential Geometry in Statistical Interference, IMS-Lecture Notes 
10 (1987) 217. 

[51] J. Lin, IEEE Trans. Inf. Theory 371 (1991) 145. 

[52] A. Majtey, P. W. Lamberti, M. T. Martin, and A. Plastino, Eur. Phys. J. D 
32 (2005) 413. 

[53] S. E. Massen, V. P. Psonis, and A. N. Antonov Int. J. Mod Phys. E 14 (2005) 
1251. 

[54] M. Edwards, K. Burnett, Phys. Rev. A 51 (1995) 1382. 

[55] E. Cerboneschi, R. Manneha, E. Arimondo, and L. Salasnich, Phys. Lett. A, 
249 (1998) 495. 

[56] P. Schuck, and X. Vinas, Phys. Rev. A 61 (2000) 043603. 

[57] S. R. Gadrc, and R. D. Bendale, Phys. Rev. A 36 (1987) 1932. 

[58] S. R. Cadre, S. B. Sears, S. J. Chacravorty, and R. D. Bendale, Phys. Rev. A 
32 (1985) 2602. 

[59] S. E. Massen, and C. P. Panos, Phys. Lett. A 280 (2001) 65. 

39 



[60] C. P. Panos, Phys. Lett. A 289 (2001) 287. 

[61] P. T. Landsberg, and P. T. Shiner Phys. Lett. A 245 (1998) 228. 



40 



